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jp 1-0 Introduction 

f ' 

pp This research proposes to use model-based control methods in the charging and discharging of 

g$ batteries and fuel cells to provide improved performance and longer life. 

m ' 

Data from previous studies of NiH2 batteries is used to facilitate the research. Based on this data, a 
low order hybrid battery model is developed This model contains some fractional order components due to 
the approximation of many diffusion processes in the battery. This model is then useable as a real-time 
observer. Damage mechanisms are studied for NiH2 batteries, which leads to the development of an 
instantaneous damage model. Using this knowledge, a control system is designed to extend battery life. 

Fundamental to the success of this research is the development of an accurate and versatile dynamic 
model of the specific battery type. This is discussed in the next section. Damage modeling is then considered 
for the purpose of developing a control system that will extend battery lifetime. Advanced versions of a real- 
time observer and observer-based control, which include features such as model-parameter tracking, specific 
battery "learning," and damage mitigation, are discussed. 

The experience gained with the MH2 batteries should be extendable to Li based batteries and fuel 
cells. Similar procedures can be followed. It is expected that application of control methods to fhel cells 
should provide some immediate durability improvements, while the life of Li based batteries can be extended. 

This research is highly interdisciplinary, as it combines the areas of electrochemistry, space power, 
control systems, dynamic modeling, fractional calculus, damage modeling, life-extending control, and power 
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electronics. An important consequence of this research is the ongoing information exchange between (he 
various disciplines involved Interdisciplinary research, such as ibis, which combines control systems and 
electrochemistry, is always fruitful, as it provides all parties with new and important insights into their fields 
of interest It is believed that improved batteiy performance and extended lifetime can be obtained through 
operational improvements by the mechanism of optimal control design. 



2.1 NiH2 Cells 

Nickel-hydrogen cells have been used in space applications for many years. The main advantages of 
nickel-hydrogen cells are high specific energy, long cycle life, long functional lifetime, and feirly robust 
chemistry. The long cycle life has been used to advantage in low earth orbit (LEO) applications where a 
satellite will orbit the earth every 90 minutes. Roughly 35 minutes of this orbit is spent in the earth's shadow, 
during which time the satellite must use the batteries for power. The other 55 minutes of the cycle is spent in 
sunlight, during which time the satellite is powered by solar panels that must also recharge the batteries. In 
these applications, up to 50000 charge-discharge cycles are common for nickel-hydrogen cells. See Linden 
for more discussion of these cells. 

Nickel-hydrogen cells are essentially a hybrid of standard cells and fuel cells. The positive electrode 
is nickel (NiOOH), and the negative electrode is hydrogen that diffuses into a porous platinum electrode. 
The platinum only serves to provide a reaction site for the hydrogen gas, and does not participate in the 
reaction. The entire cell consists of a large sealed pressure vessel that contains the hydrogen gas and the 
electrode stack. Hie electrode stack alternates pineapple-slice shaped disks of nickel electrode, separator 
(containing the electrolyte, KOH), platinum electrode, and a hydrogen diffusion screen. Collections of cells 
form abattery as shown in Figure 2.1. 



2.0 NiH2 Batteries 




Figure 2.1 



Typical NiH2 battery configuration (Briton & Miller). 
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2.2 NiH2 Cell Chemistry 

In normal charge-discharge operation of a nickel-hydrogen cell, the electrochemical reactions are 
as follows: 



(2.1) 



negative hydrogen electrode: £H 2 + OH~ H 2 0 + e~ . (2.2) 



positive nickel electrode: NiOOH+H 2 0 + e~ Ni(OH)2 + OH 



Durmg discharge at the positive electrode, the NiOOH (nickel oxy-hydroxide) on the electrode reacts with 
water in the electrolyte at the surfece of the electrode as it accepts an electron from the external circuit 
through the nickel substrate. The products of this reaction are the Ni(OH)2 (nickel hydroxide) on the 
electrode, and OH" (hydroxyl ion) which diffuses away from the from the electrode through the electrolyte in 
the separator. At the negative electrode, the hydrogen gas reacts with hydroxyl ions in the electrolyte to form 
water and release an electron in to the external circuit This process is reversed during charging. The 
electrolyte, which is often 31% KOH solution in water, contains potassium ions, K + , and hydroxyl ions, OH", 
dissolved in water. Clearly the water and the hydroxyl ions are important in the reactions at the electrodes. It 
is usually assumed that the potassium ions do not take an active part in the cell reactions, however its exact 
role in the overall cell behavior is unclear. 

During char ging, overcharge reactions can occur. These occur whenever the charging current drives 
the above reactions fester than the required mass transfer can happen. This can occur both when the nickel 
electrode is fully charged, and when the charging current is too large. Thus, in overcharge, the nickel and the 
hydrogen do not participate in the reactions, and the electrolyte is reacting electrochemically at both 
electrodes. The overcharge reactions are as follows: 

positive nickel electrode: 20H" ove ^° / * g > H 2 0+£0 2 + 2e~ (2.3) 

negative hydrogen electrode: H 2 0 + £0 2 + 2e~ — overcharge > 2 OH" . (2.4) 

Consideration of these reactions would seem to indicate that overcharge is not a problem, but the formation of 
oxygen in a hydrogen environment is potentially very dangerous. This is discussed at greater length in the 
damage section below. 

In a collection of cells composing battery, cell reversal can sometimes occur when a single cell 
discharges completely before other cells. This cell then becomes a load on the other cells. If the cell is 
assumed to be nickel limited, the reversal reactions are as follows: 



positive nickel electrode: H 2 0 + e" rgvg ™ / > OH" + -| H 2 (2.5) 

negative hydrogen electrode: £ H 2 + OHT reversal > H 2 0 + e " . (2,6) 

It should be observed that the hydrogen gas that is generated at the positive electrode is consumed at the same 
rate at the negative electrode. Thus there is no net increase in hydrogen pressure or any net change in 
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electrolyte concentrations during reversal, and the cell can be operated continuously in this mode Cell 
reversal is not considered further in this paper. * 

3.0 Performance Model for NiH2 Batteries 

3.1 Dynamic Modeling Background 

Control systems are often used to improve the performance of many systems and devices (Friedland 
(1985)). Typically die quality of the control is related to the level of understanding of the system being 
controUed(Dorf(1974)). This system knowledge is contained the system dynamic model. Useful models for 
control design are dynamic, that is represented by differential equations, and are of reasonably low order 
(requiring only a few differential equations) while still capturing the important dynamics of the specific 
system (Hartley, Veillette, DeAbreu, Chicatelli, and Hartmann (1998)). A more complicated model, such as a 
computational fluid dynamic (CFD) model, will have thousands of differential equations (or dynamic states). 
This type of model is useful for a detailed understanding of battery phenomena and for cell design. Such a 
model, however, can contain states that are impractical to control, and can lead to a control system that is 
m unnecessarily complex (Chicatelli and Hartley (1998)). Furthermore, a CFD model is often only represented 
R *s a computer code, and effectively no equations are available to the control system designer. On the other 
W hand, a less complicated model, such as an empirical input-output relationship, may not contain any 
m differential equations at all. This static model will usually not allow good control system design. Thus, to 
m design a good control system, a low order dynamic model of the system to be controlled is most useful This 
up. model should capture enough of the dynamic behavior of the system to be accurate over the frequency range 
jsll of interest, while not being overly complicated. 

ib In the area of battery and fuel cell modeling, there would appear to be two common types of model, 

Q CFD-based and empirical. The highly accurate CFD-based models are based on first principles, and 
0 accurately model the diffusion of all of the active species through whatever media are present (Doyle, Fuller, 
£3 and Newman(1993), Newman(1991)). Such models require a large number of electrochemical parameters, 
^ including things such as ion mobilities and reaction rate constants. Also required is specific geometrical 
Q construction data on a given cell. These models are highly accurate and provide good representations of cell 
fif dynamics as long as the cell configuration remains constant Unfortunately, as actual cells are charged and 
discharged, their dynamic parameters do not remain constant, and these models typically do not provide good 
representations for electrochemical cells over many charge-discharge cycles. Additionally, the exact 
construction varies from cell to cell, and thus these models do not quite have the actual cell information. In 
terms of using these models for control, there are also a large number of states due to the representation of 
many species diffusing through the cell. Furthermore, the spatial lumping of such diffusive systems usually 
yields a dynamic system that is very stiff (Hartley, Beale, and Chicatelli (1994)). Such systems are difficult 
to simulate, and rather than programming such a set of differential equations, an implicit simulation method is 
implemented directly on the diffusion equations in the computer code without the generation of any useful 
mathematical model on paper. Finally, these CFD-based models are expensive to build and change, while 
also being slow to build and use. 

At the other modeling extreme from the CFD-based models are the empirical models. These model 
typically have a specific algebraic structure, and the model parameters are identified from given data 
(Shepherd(1965)). These are usually simple models -that may provide reasonable static operating point 
representations of a given cell as it charges or discharges . Other than having a built in integral to track the 
charge, these models are usually not dynamic. Consequently, they can be considered as static maps, such as 
that of a transistor or of a compressor. Although they can predict an operating point, they do not accurately 
represent dynamic changes. Thus, the empirical models are not adequate for control system design. 
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« nerfn™,^! c^on ^proach to battery and fuel cell modeling is what could be termed a hybrid model, or 
ISSET^ 1 i?*? of *t s ^ to *e important first principles dynamics, while not 

be«nmng so ^heated as to not be readily useful. Their parameters are then tuned to match experimental 

^nn« Ca 5So^f e ^^ Zmg model " Such » model for electrochemical cells is diseased in Xia 
and Hartley (2000a, 2000b) This model has a chosen structure based on an electrochemical understanding of 
cell behavior. Given tins dynamic model structure, the parameters of the model are detennined for a given 
cell trom charge-discharge data. An important aspect of this model is the representation of the many 
^fusion processes internal to the cell by a fractional order system (Lorenzo and Hartley(1998)). 
Consequently the model is of low order (few dynamic states with some fractional) and can be used for many 
duTerent cell chemistries. A model of mis type is fast, understandable, and amenable to control design. 
Furthermore, it is possible to update the parameters in the model as the cell changes through its charge- 
discharge cycles. Thus the hybrid modeling approach should provide a useful model structure for 
understanding cell changes through time, while also providing a useful dynamic model for control design. 

3.2 Overview of NiH2 Cell Dynamics 

As in most electrochemical systems, the dynamics of NiH2 cells are quite complex. This section 
presents an overview of the electrical and material behavior in a cell. Figure 3.1 gives a generic structure of 
Qie major components of the cell dynamics in normal operation. Here the electrical and material contributions 
have been separated for clarity. At each electrode, a voltage is generated between the electrode and the 
electrolyte based on the concentrations of reactants. The reactants then diffuse into the reaction site or away 
from the reaction site. The hydrogen gas must diffuse into the platinum electrode from the ambient chamber. 
The nickel compounds may not diffuse, but it is not clear if the reaction front at the nickel electrode is diffuse 
or not Also, the hydroxyl ions and me water must diffuse from one electrode to the other through the 
electrolyte. The electrolyte diffusion represents a major resistance to the flow of charge through the cell. An 
additional source of electrical resistance is the charge transfer resistance at die electrode surface, across which 
there is also a double later capacitance. A voltage profile inside a typical cell is shown in Figure 3.2. There it 
is shown that voltage is produced between the electrode-electrolyte interfaces. It is shown later in this section 
that these voltages are dependent upon species concentrations on either side and on the current passing 
through Additionally there is always an ohmic voltage loss in the electrolyte due to current passing through 
It should be noted that this ohmic loss always corresponds to heating. 
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Figure 3.1 Overall NiH2 cell dynamic components. 
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Typical voltage profiles inside a cell 



3.3 Model Essentialization 

To extract the essence of the dynamics from the complex model given in Figure 3.1, several 
assumptions are made. The major assumption is that one electrode dominates the dynamics of the cell In 
this case, we are assuming a hydrogen rich environment, and the nickel electrode is dominant in the dynamic 
behavior. At this electrode, we will only keep track of the amount of charged material, c„ not the discharged 
material on the electrode. We will also lump all of the diffusive behavior into one fractional order term, c*?, 
which alleviates the need to model several diffusion channels in the cell. Finally, the voltage generated at the 
electrode will depend on the concentrations of both the stored material and the diffusing material at the 
electrode surface, and will be used as the terminal voltage. With these assumptions, the essential model is 
shown in Figure 33. 

In the Figure 3.3, the electrical behavior is shown in the rigftt side of the figure, and the chemical 
species behavior is shown on the left side. On the electrical side, the voltage at the terminals, v, is due to the 
dependent voltage source that is a function of electrical current, i, concentration of charged material, c S9 and 
concentration of diffusing material, at the electrode. On the chemical side, we are using an electrical 
analogy of the material mass transfer. In the analogy, mass flow rate becomes current (or flux), j, and 
material concentration becomes voltage. Thus the two dependent current sources on the chemical side are 
mass flow rates, j a andj d , defined by the current on electrical side. The charged material is stored in a 
capacitor, C, that represents mass storage, with the addition of the self-discharge resistor, JRsd, representing an 
orifice for the loss of charged material. The diffusing material must pass through a fractional-order element, 
fo"* that is an agglomeration of all the various diffusion pathways in the actual cell This fractional-order 
element separates the concentration of diffusing material at the electrode, c d , from the concentration in the 
bulk of the electrolyte, c^ 
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The effective chemical reaction in the cell then becomes 



S+D + e" reaction products , (3.1) 



charge 



where charge conservation is assumed in the essentialized model. Here S represents the reacting stored 
material and D represents the reacting diffusing material. We will often associate the stored material in a 
NiH2 cell with NiOOH, and the diffusing material with either HjO or OH", but this is only for mental imagery 
and not inherent to the math model We could also keep track of the reaction products if it were necessary for 
the accuracy of the model, but we have not needed to do that in our study. Based on this approximate cell 
reaction, during discharge, the stored material reacts with the diffusing material and an electron. Thus when 
current is leaving the battery (electrons entering), both the stored and diffusing material are being deleted at 
the active electrode. Note that we could have assumed the diffusing material was building up outside the 
electrode during discharge, rather than going to zero, but that would have added a little more complexity. 

. The next several sections discuss the specifics of the various elements of the model. A discussion of 
toe idennfication of the specific constants in the model is included. Note that a thermal model is not 
developed directly for the essentialized model. It is assumed mat the heating effects are related to the cell 
damage, and that the damage model developed in Section 4.4 will contain the necessary information. Heating 
is considered further m Section 4.3 where energy flow, or power, in the cell is discussed 
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3.4 Electrode Behavior: Faraday's Law 

For the discussion of the next two sections, consider the electrode process : 

oA+6B< > cC+dD+ne~. (3,2) 

This states that chemical species A and B are converted to species C and D plus electrons, at rates related to 
the coefficienls a,b,c,d,n (it is also assumed that these compounds are appropriately ionized so that charge is 
conserved). It is important to understand that this equation represents both a chemical process and an 
electrical process. In feet, the reaction rates can be completely determined using the electrical process by 
recognizing that the chemical conversion can only occur if electrons are either arriving or leaving. Thus, the 
chemical conversion rates are controlled by, or measured by, the electrical current passing through a given 
electrode. All the rates of chemical conversion are then known via the following observation. For every n 
electrons leaving die electrode, c molecules of C are produced, d molecules of D are produced, a molecules of 
A are consumed, and b molecules of B are consumed Then recognizing that the rate of electron production is 
related to electrical current /, the following rate equations result : 

* n dt n dt n dt n dt 9 ( } 

Where q r is the charge of a single electron. It should be noted that concentrations of chemical species are 
often used Concentration is defined as moles of the species per unit volume. Concentrations are also often 
represented by [A] or by ca, which would represent the concentration of species A. Using concentrations, the 
previous equation becomes 

.i = „F% = .A = «A^ + a^ = + ^ (3.4) 
dt dt dt dt dt V ' 

where Fis the Faraday constant relating the moles of electrons in a coulomb, F=96484.6 Coul/mole. 

3.5 Electrode Behavior: Electrode Equation 

Once the concentration rates on either side of the electrode interface are known, it is important to 
obtain an expression of the voltage in terms of these concentrations. This can be done by considering 
conservation of charge through the interface. Again assuming the reaction of Equation (3.2), the required 
relationship is determined by considering the fluxes, J 9 of the species in the forward and reverse reactions at 
the electrode, 

*=Jf~Jr- (3.5) 

Assuming that the species fluxes are proportional to concentrations, yields 

i = k f c* A c b B -k r c c c c%. (3.6) 

Here, the rate constants, A, can be related to the electrical potential across the electrode-electrolyte interface 
using free energy considerations, 
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k f = k 0 e iFlBT) ^^ (3.7) 

K-ko^'*™* 9 *^ (3.8) 

where R is the universal gas constant, Tis the temperature in K, or is an electrode energy barrier symmetry 
fector, v - v 0 is the voltage deviation from some nominal voltage v 0 , and ko is a rate fector. 

Inserting these gives the electrode equation, 

/ = * 0 {c°c> e <r - c c c c% e<™^) ) . (3 9) 

This equation is important to both the electrical and chemical behavior of the electrode-electrolyte interface as 
it relates the voltage across the interface to the concentrations of active species on either side of the interface 
and to the electrical current passing through the interface. It is often called the current-potential 
characteristic, and the approach used is often referred to as then Butler-Volmer approach (Bard and Faulkner 
(1980)). 

Thus, for an electrode-electrolyte interface satisfying a reaction given by Equation (3.2), the required 
equations relating electrical current through the interface, electric potential across the interface, and species 
concentrations on either side of the interface are Equation (3.4) and Equation (3.9). These equations can also 
be used in a multiport configuration to relate and separate the electrical variables from the chemical variables. 
Consequently, these two equations are fundamental to representing the behavior of any electrified interface. 

3.6 Linear-in-the-parameters Electrode Equation 

The electrode equation is highly nonlinear, particularly when solving for the voltage in terms of the 
concentrations and the current. It is almost always solved using numerical methods, such as Newton — 
Raphson iteration. As an alternative to doing this, we propose a linear-in-the-parameters approximate 
solution to the electrode equation. The form we have chosen is 

v=fc t +Jc 2 ln(l+| z|)sgn(0+A 3 ln(c rf )+* 4 ln(l -c,) (3.10) 

where the variables i 9 c s ,c d correspond to those in our essential model and the &'s are parameters to be 
determined from data. What this equation assumes is that the terminal voltage v, is equal to the sum of four 
terms, and is effectively a compromise between the Tafel and the Nernst solutions of the electrode equation 
(Bard and Faulkner). The Nernst solution assumes that the current is so small as to be negligible. The 
electrode equation is then easily solved for the voltage in terms of the concentrations, 



v = v 0+ gln| - (3.11) 



Alternatively, the Tafel solution assumes that the current is large in one direction or the other, which means 
that one of the two exponential terms is negligible. With one of the exponentials equal to zero, the electrode 
equation is easily solved for the terminal voltage in terms of the current, 

V = V ° + " J^F (3.12) 
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In the linear-in-the-parameters approximate solution to the electrode equation, Equation (3.10), the 
first term is the constant which represents the voltage at the fully charged, zero current, steady state 
concentrations condition. It also represents the v 0 terms in the Nernst and Tafel solutions. The second term 
represents the nonlinear logarithmic dependence on current from the Tafel solution. Hie problem with the 
Tafel solution is that it assumes the current is large in one direction, but not both, while our approximation 
must be applicable to currents in both directions. Taking this into consideration yields our compromise term 
with the abs and sgn functions. The third and fourth terms represent the Nernst dependence on concentrations 
. on either side of the charge storage reaction (and ignoring the feet that one reactant is in fee solid state). The 
third term represents all the collective diffusion tarns, as well as the electrolyte ions that are being consumed 
or generated at the electrode surfece. As this concentration goes to zero, the terminal voltage gets small. As 
this concentration goes up, it has less and less effect on the voltage. The fourth and last term represents the 
effect of stored charge. As the stored charge approaches one, the voltage gets big. If it ever gets to zero, the 
voltage will also get small. 

Another major benefit to our linear-in-the-parameters approximate solution is that it is possible to 
■ identify the 4 ^-parameters directly from charge-discharge data using the batch or recursive least squares 
method (Goodwin and Sin). The requirement to do this is the actual measured terminal voltage and current, 
along with the estimated transient values of the stored charge and the diffusive material 

Note also that the charge transfer resistance through the double layer is contained implicitly in the 
current term in the linear-in-the-parameters approximation. Likewise, we are not modeling the double-layer 
capacitance explicitly, as it is also implicit in the linear-in-the-parameters approximation as well as the 
diffusive term. 



3.7 Charge Storage 

The charge storage will be through charge accumulation by integrating the current Using the 
electrical analogy, this element is effectively capacitive. Any charge that passes into or out of the cell, will 
directly effect the concentration of the stored charge via Faraday's law, 

4,0) -W - (3.13) 



3.8 Diffusion Behavior 

The diffusive components are lumped into one fractional-order term. This term effectively represents 
the diffusion of material toward or away fft>m the electrode through the electrolyte, and with our assumed 
reaction, is thus on the same side of the reaction equation as the stored charge 



S + D + e~ ^ reaction products. 



As a fractional-order element, its transfer function is 
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where s is the Laplace variable, and kj and q are constant to be determined from data. However, the diffusive 
material concentration must return to the assumed bulk value when the system is quiescent Thus die 
equation is converted to the time-domain and modified to reflect this return to bulk concentration 

** CO =<W - k d 0 d;«i(t) . (3. 15) 

It should be remembered that fractional-order terms have time-varying initialization functions associated with 
them. These were determined to have an overall effect on the diffusive terms of less than 3% our specific 
study, and were not included (ignored). If the implemented battery model is used in other ways after this 
study, the initialization functions may need to be included, which might also require further adjustment of the 
constants in Equation (3.15). 

3.9 Self-discharge 

Self-discharge is implemented through a simple resistive (to flow) loss across the charge storage 
element, 



f «* - 4*- • (3.16) 

K sd 



This should be thou^rt of as a pathway for stored material to leak out of storage (bigger less loss). 

3.10 Ohmic Cell Resistance 

Most conducting media provide some resistance to the passage of electrical current. We could ' 
attempt to use one of the many ways that this is often calculated in electrochemical cells for our model, but 
we will allow the current term in the linear-in-the-parameters electrode equation to implicitly represent this 
resistance. If a value of resistance is needed, it can be obtained by linearizing the linear-in-the-parameters 
electrode equation. We will do that here since the resistance of a cell is often talked discussed. Noting that 
resistance can be defined from Ohm's law as 

n dv 

lu' CUT) 

then 

R _d V _ d(k 1 +kM^+U\)^i)+k 3 bi(c d )+k 4 ]n^-c s )) 

Re*** - -ji Ji ' < 318 > 

or 

dv _ <*(* 2 ln(Hf|) Sgn (Q) 

K ohmic ' (3.19) 

We can evaluate this for any steady-state current Evaluating at zero steady-state current for symmetry gives 




(3.20) 
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^n^r™ e "f a Tf IS ? le f trode equation &™ a measure of the total Ohmic resistance in 
^Sl^ b0tb fte Ch3rge 81 tetrodes, and the Ohxnic 



3.11 Essentialized Model Overview 



highUgh^e^uto^L 8 ^ * KSUre ^ ™ e SySt6m *— * eqUati0nS ^ ***** beIow to 

Tenninalbehavion current into the battery Is +f, (3.21) 

the terminal voltage is +v, ' 



stored material with self-discharge: c, (r) = i(r) - — c s (/) , (3.22) 

diffiising material: ^(0=^-^0^/(0, (3.23) 

electrode equation: v(c„c d ,0=k 1 +k 2 mfc + |f|)*gn(0+*s ln(c rf )+A 4 m(l-c,) . (324) 

From this we can see that the required parameters are, self-discharge resistance R, d , fractional-order q and 
constant k* bulk concentration of the diffusing material c bulk , and the linear-in-the-parameters electrode 
equation constants, A,, fe, fe, and A4. 



3.12 Parameter Determination 

,. ^. To most eas * r discuss model parameter determination approach, we will apply it to a specific 
cell. The cell chosen is the NSWC Crane Pack ED 3602G (Gates). It is rated at 65 AHr, and uses 31% KOH 
concentration, and is maintained at 10 degrees C. The charge-discharge profile is a square wave, with 35% 
depth-of-discharge (DOD) and 104% recharge ratio. The square wave current is 26.29 A for 54 minutes 
charging, and -37.92 A for 36 minutes during discharge. Note that 65 AHr = 3900 AMin, 35% of 65 AHr = 
1365.1 AMin = 26.29A*54Min, 3900AMin - 1365.1 AMin = 2534.9 AMin, and 104% of 1365.1 AMin 
=1419.7 AMin = 37.92A*36Min. 



The self-discharge resistance is determined by the recharge ratio of 104%. We will assume that 
roughly half of the extra 4% goes into self-discharge, and that the other 2% goes into overcharge reactions. 
The division of the excess could be determined more accurately using the Buder method (Buder, 1972), but 
we have not done this at this time. Assuming no current in or out, we assume mat 2% of the stored material 
(1.3 AHr = 78 AMin) is lost during each 90 Min cycle via self discharge. Using this in the known solution to 
Equation (3.22) 

«•<©-«, COD*"*'*' 

gives 

c, (0 =3900 - 78 = 3900e _90/ **' . 



(3.25) 
(3.26) 
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Solving this for the self-discharge resistance gives R^ =4455 ohms. This value was later modified after the 
simulation was implemented to force the constant-plus-taper charging to give a final charge of 3900AMul 
The modified value of the self-discharge resistance that is used in the charging studies is R sd =4024 ohms, or 
equivalents 1/R sd =0,0002485 . 

To determine q and kj for the diffusive material, we will observe that most NiH2 cells are diffusion 
limited. This must be true, since the terminal voltages get small after only 35% of the active stored material is 
used. To start, it is assumed that the bulk concentration is unity, 



^buBc 



= 1. (3.27) 



Thus we will assume that the diffusive material goes from unity concentration to roughly zero concentration 
as 35% of the stored material is used. With a constant current of -37.92 A, the solution to the diffusing 
material equation using Laplace transforms is 

The time response obtained by inverse transforming this is 

c rf (0=l- 37.92V*. (3.29) 

To determine the unknown constants, we will evaluate this at the end of discharge, or *=36 Min, at which 
point the concentration goes to roughly zero, 

1-37.92*, (36)* =0, 

or 

37.92*, (36)* =1. (3.30) 

Taking the natural logarithm of both sides and rearranging gives 

-ln(37.92) = gln(36)+ln(Jfc rf ) . (3.31) 

This gives one equation for our two unknowns. If we now assume that the charging occurs similarly, we can 
obtain another equation. This is not exactly true, since the diffusive material will recover on its own, with no 
driving current, but this recovery is somewhat slower than the forced recovery. We will thus ignore the 
associated initialization response during the charging phase. Thus, charging at 26.29 A for 54 Min yields the 
equation 

- ln(26.29) = q ln(54)+ln(*„ ) . (3.32) 
These two equations can now be solved for the two unknowns, which become ^=0.9034 and fe=0.001036. 

To implement the diffusion term, an approximation to the fractional operator is necessary. We will 
use the approximation of Sun, Abdelwahab, and Onaral (1984). The resulting approximation is 

A„-~ = 0.001036 0-01f« 2 +l-1444'+(>-345S 
" s + 0.5078s +0.0013 
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which is designed with a maximum Bode magnitude response error of 2dB between m = lO 3 and a> - 10 +2 . 
Better approximations containing more terms, or simpler approximations containing fewer terms are easily 
generated using the method. Again it should be noted that the initialization terms were ignored in this study. 

The parameters in the electrode equation are obtained by using a batch least squares optimization on 
the charge-discharge data from the actual cell (see program in appendix). The resulting equation is 

v=13656+ 0.026Sln(U|i|^^ (3 34 } 

Notice that the last term has been rescaled so feat c 3 can go to the AHr capacity of the battery which is 
3900AMin. ~ 7 



3.12 Model Verification Studies 

Figures 3.4, 3.5, and 3.6 show actual charge-discharge data (the 'o* points) from the cell discussed in 
Section 3.12. These figures also contain the charge-discharge plot (solid line) from the identified model also 
discussed in Section 3.12. The three plots are chosen at 500 cycles, 10000 cycles, and 40000 cycles 
respectively to shown how the specific cell has changed with time, and to show the ability of the model to 
track changes in the cell. 




Figure 3.4. Charge-discharge cycle for the chosen cell at the 500 th cycle. 



5/7/02 9:08 AM 
14 



Optimal Battery Charging for Damage Mitigation 



Tom T. Hartley and Carl F. Lorenzo 




Figures 3.5 and 3.6. Charge-discharge cycle for the chosen cell at the 10000 and 40000 cycle. 
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4.0 Damage Modeling in NiH2 Batteries 

4.1 Background 

If the real-time observer is going to be able to provide information concerning the state of cell health, 
some understanding of the damage process in an electrochemical cell is necessary. There is no known 
mstotaneousj damage rate model for electrochemical systems, however, fetigue and creep damage modeling 
has been used for damage mitigating control in aerospace propulsion systems (Lorenzo (1994)). In that 
application, damage per cycle information has been used, as well as a continuum damage model developed 
The continuum damage model predicts the current damage rate as a function of local stress and strain, and can 
be used to implement a controller that trades a small amount of performance (dynamic response) for large 
gains in system lifetime. A continuum damage model will be required for work with batteries and fuel cells. 

For electrochemical cells and batteries, a number of damage mechanisms have been identified which 
include: corrosion, electrolyte communication between cells, separator dry out, popping damage, plaque 
expansion and electrochemical degradation, Still other mechanisms may be present which have not yet been 
identified. However it still remains to quantify these mechanisms, determine their importance to a particular 
design, and to determine the operational relationships to damage accumulation. These activities are beyond 
the scope of the current research, however as these effects are quantified, the resulting relationships can be 
used to refine battery damage prediction. The next two sections discuss in a general way the mechanisms that 
may be responsible for battery cyclic damage. 

4.2 Damage Mechanisms for NiH2 Batteries 

There are many ways to damage any cell. We will focus on the damage mechanisms that we think 
are presently the most significant for NiH2 cells. 

Formation of y - phase NiOOH: 

At the positive nickel electrode, the charged material is NiOOH, and the discharged material is 
Ni(OH) 2 . As in many solids, however, these compounds can exists in difference solid phases. The normal 
charge-discharge cycling occurs between the beta phases, 



/?-NiOOH /?-Ni(OH)2 • 

charge 

On overcharge, however, the excess heating due to O2 formation along with any large currents, causes the fi - 
NiOOH to undergo a phase change into the y -phase, y -NiOOH possesses two properties of interest to our 
studies. One is that it generates a lower voltage with respect to the KOH electrolyte, roughly IV and opposed 
to 1.5V for fi -NiOOH. The other physical property is that the spaces between the molecules are much larger 
for y -NiOOH than for fi -NiOOH. Thus the formation of y -NiOOH can physically damage a cell stack due 
to the changing size associated with the phase change. Additionally, as the amount of y -NiOOH builds up, 
the usable charge storage capacity of the cell goes down, as the y -NiOOH cannot provide the voltage 
necessary for the electrical loads. 
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Cells with a build-up of y -NiOOH can be reconditioned using the reaction mechanism proposed by 
Bode. This is shown in Figure 4. 1. It should be observed that the y -NiOOH will cycle electrically with a - 
Ni(OH>2. We can recover the original fi -Ni(OH) 2 by allowing the cell to stand for an extended period of 
time in a KOH solution. 



fi -NiOOH 



overcharge 



-> y -NiOOH 



charge 



charge 



discharge 



/?-Ni(OH) 2 



a-Ni(OH>2 



standing in KOH 
Figure 4.1 Solid phase relationships for a NiOOH electrode. 



Formation of 02: 

Overcharging a cell can occur by either continuing to charge the cell after all the fi -Ni(OH) 2 has 
been converted to fi -NiOOH, or when die charging current is so large that the fi -Ni(OH>2 cannot react fast 
enough to form fi -NiOOH. In either case, the excess current goes into the electrolysis of the cell electrolyte. 
The effect of this is the formation of 0 2 at the nickel electrode, along with heating. The problem is further 
complicated by die formation of oxygen in a hydrogen rich environment. These can recombine chemically 
and release significant heat This recombination can also do physical damage to the electrode. This chemical 
recombination of oxygen and hydrogen is often called popping. Clearly, overcharging is to be avoided. 

Heating: 

Heating is a major problem in NiH2 batteries. As discussed above, it results in the formation of y - 
NiOOH, which reduces cell capacity and does physical damage to the cell. Another problem with heating is 
the increase in self-discharge reaction rates, which can also be manifested as a reduction in cell capacity. Any 
excess heat production must be avoided. 



43 Thermal Analysis and Hysteresis 

This section will discuss the gross energy balances taking place in an electrochemical reaction. We 
will start with the free energy of the reaction, 

AG = AH - T6S . ( 41 ) 
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Here AG is the Gibbs free energy of the reaction and represents the energy available for external work, and 
in electrochemical cells is effectively the source of the terminal voltage, Aff is the enthalpy of the reaction 
and represents to total chemical energy change during the reaction. The TAS term is the temperature 
multiplied by the entropy change of the reaction. The entropy change involves the reordering of the reacting 
molecules as they undergo chemical reaction. Multiplying the entropy change by the temperature of the 
reactants, TAS , gives the heat released or absorbed during the reaction. If the free energy is negative, the 
reaction will be spontaneous. 

Dividing the above equation by -nF, where n is the number of electrons involved in the reaction and F 
is Faraday's constant, gives 



AG _ AH j TAS 
nF nF nF 



+ -^- (4.2) 



The sign convention of a spontaneous reaction giving a positive voltage in an electrochemical cell then gives 

AG = -nFv. (43) 
Using this in the energy equation yields 



TAS 
nF 9 



Here 

AH = -nFv m q (4.5) 

is the thermoneutral voltage, or the voltage of the electrode reaction when there is no heat being generated. 
Although free energies, entropies, and enthalpies are all generally dependent on the concentrations of 
reactants, the enthalpy of reaction for many electrochemical reactions is roughly constant This allows the 
thermoneutral voltage to be used as a measure of when there is heating at the reaction site by comparing it 
with the actual terminal voltage. 

The energy flow in the reaction can be obtained by multiplying the above equation by the current, 

or 

rv = rv w +e r - (4.7) 

Here Q r is the rate of heat release or absorption at the reaction site, and represents the unavailable electrical 
power. This equation can be rearranged to give 

Qr = '"(v-v^) (4.8) 

which shows that comparing the terminal voltage with the thermoneutral voltage allows determination of the 
heating due to the electrochemical reaction. 
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The above analysis has only considered the energy related to the electrochemical reaction, and has not 
considered any of the electrical behavior elsewhere in the overall cell, particularly any electrical losses. 
Whenever current flows in a conductor, there are electrical losses due to the heating of the material. We will 
use the total Ohmic resistance as the primary electrical heating source in this discussion. Doing so, we can 
write Equation (4.4) as 

v = v w +^ + « ofani ,. (4.9) 

This shows that the terminal voltage varies not only with the free energy of the reaction, but also with the 
. resistive loss in the cell. Multiplying this equation by the current gives 

iv = iv m +i^ + i 2 R olmlc . (4.10) 

This equation can be expressed in words as: 
01 electrical power = electrical power of the reaction + reaction heatingteooling + electrical heating, 

fjjj Thus the overall heating is a combination of electrical and chemical effects. 

m 

Hf . . The energy flow can be further generalized as in Sanabria, by including the energy that goes into 

§|f oxygen production during overcharge. Equation (4. 1 0) is generalized adding this electrolyzing current, 

7 iv = {i-i m > m Hi-ioz)™ + ioziv-iRoHiae) + ***** • ( 411 > 

Q 

HI In words, this says: 



electrical power in = electrical power into the reaction + reaction heating/cooling 

4- heating due to 02 production + electrical heating. 

The discussion of this section can be summarized in Figures 4.2 and 4.3 which show charge-discharge 
curves onav^ diagram, along with the associated energy distribution as described by Equation (4. 10). These 
diagrams are somewhat ideal in the sense that they do not have self-discharge or overcharge effects included 
In the case of overcharge, the figure would need to be modified to reflect the current going to the oxygen 
production as described in Equation 4.11. Another important effect that is not considered in these figures, or 
the discussion above, is that the voltage can move up or down as the diffusing material relaxes (independent 
of the stored charge). Since only the stored charge is plotted on the horizontal axis, the figures do not take the 
concentration of diffusing material into effect 

At this point, sufficient data does not exist to quantify these effects for a particular battery design. 
Thus, based on available cycle life vs depth of discharge data we have developed an empirical continuum 
damage model for cells based on external variables, such as voltage and current, and internal variable such as 
species concentrations. It may be possible to explicitly include other internal variables, such as pressure and 
temperature in the future. This damage model allows the development of a damage mitigating control for 
NiH2 cells which is discussed in the next sectioa 
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Figure 4.2 A typical charge-discharge 
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cycle, showing the energy distribution during charging. 
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4.4 Battery Continuum Damage Modeling 

This section develops methods for using battery cyclic life data to infer instantaneous damage rate 
within the battery. This ins ta nt a ne ous damage rate is required to design optimum charging profiles and to 
establish control algorithms. Cycle life as a function of depth of discharge (DOD) data for NIH2 cells is 
presented in the Green-Hoffman report The following data from that report is given for a temperature of 10 
degrees C, 



DOD 


Cycles 

To 
Failure 


35% 


38,000 


50% 


19,000 



For this data, a charge/discharge cycle has a 90-minute period, with a 54 minute charging duration. The 
recharge ratio is always less than 1 10%, and the recharge current was constant until 94% recharge at which 
point a taper charge was used to finish the charging. The cycle life data is reproduced in Figure 4.3. In this 
figure, only two data points are available at 10 degrees C, Four curves have been fit to the data, the Green- 
Hoffman model, a modified Green-Hofiman model discussed later, a power law, and a linear curve. 
Although only the two Green-Hofiman approaches are considered, airy curve can be fit through two points. 
The Green-Hoffman approach uses the following exponential model 

N fcH «1885.04« WM<WIQD) , atl0°C (4.12) 

where N f is the number of cycles to failure. We will assume that when the total accumulated damage 
becomes equal to one, the cell has failed. Otherwise stated, the total accumulated damage is normalized to 
one. Further assuming that the damage per cycle is constant for a given DOD, or more specifically, that the 
damage per cycle does not increase as the accumulated damage increases, then the damage per cycle is 

A** = O.OOOSdOSe- 4 * 21 =—5—. (4.13) 

f GH 

Based on similar previous studies of mechanical fatigue (Lorenzo), we will derive an expression for 
the instantaneous damage, D(t) , based on the particular cycle life model This is possible because the cycle 
life data implicitly contains all the information about how heating, overcharge, high current, and other hard to 
model factors damage a cell. The cyclic damage model expressed by equations (4.12) and (4.13) 
characterizes the damage hysteresis cycle size by DOD, which is actually representative of tire charge 
amplitude of the cycle. The hysteresis cycle size could just as well be characterized by the voltage amplitude 
of the cycle. The purpose of the models to be derived here is to optimize the charging current profile. Models 
based on DOD as characterizing the size of the damage cycle are thus not useful because the overall charge to 
be stored during charging is always the same, no matter what charging profile is used An alternative and 
more appropriate descriptor of cycle size for damage modeling is the voltage amplitude of the cycle. The 
voltage amplitude referenced to the minimum voltage point of the cycle (i.e. the discharge cutoff voltage 
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(35%)) will be used to describe the damage cycle. The stress variable then will be the voltage difference 
based on this reference voltage. 




Figure 4.3 Cycle life vs DOD data. 



For a full charge/discharge cycle, the voltage amplitude v a varies with DOD (charge amplitude). This 
relationship is certainly a nonlinear function, but we will use a linear approximation. We express this as 

DOD=c t v a . (4.14) 

To obtain c l9 we adjust it until the simulated damage per cycle for a constant-plus-taper charge is equal to 
that given by Equation (4.13) for 35% DOD. This equation is substituted back into the damage per cycle, 
equation (4.13), giving 

= 0.0005305 e- 4 ' 621 * 4 - 621 ^ . (4.15) 

Now the accumulated damage as we traverse a charge/discharge cycle is the damage per cycle. Thus the 
integration of the instantaneous damage through the cycle is 
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fav)dv = D^ = 0.0005305e-* 62,+ * 6a,e » v <' . (416) 



cycle 

where S(y) is the voltage referred damage rate, that is 

^(v) = — . (4.17) 

We now assume that all of the damage is accumulated over the charging part of the cycle, that is die "uphill" 
part ofthe process, then 

p(v)rfv = = O.OOOSBOSe^ 1 * 4 * 621 ™ = 5.222 xKTV' 621 ^ . (4.18) 

Vmin 

We note that models including damage on the discharging portion of the cycle are readily possible using 
techniques similar to those used here. Changing the limits of integration gives 

p(v + v^)dv m O.mSMSe-** 21 * 4 * 11 ™ m 5.222 xlO^e 4 -* 21 ™ . (4.19) 
o 

Since v a =v max -v^ , we can now infer an approximation to the shifted voltage referred damage rate, 
^(v + v^), as 

<?(v + v^) s OMQSZOSe** 11 4.621^™* v (420) 

The reason that this is approximate is due to the fact that at cyclic amplitude of DOD- v a =0 , the Green- 
Hofl&nan model does not yield zero cyclic-damage, as we require. In feet, as discussed in Section 4.2, zero 
amplitude would actually allow healing. Then Equation (4.20) becomes 

S(y) =2.4129x I0^c x e 4 ^^^ . (4.21) 
This result in Equation (4.18) gives 

v raax 

jS(y)dv= 5.222xl0- 6 (e 4 - 62,<!,v « - 1)= 5.222 x 10"* e 4 - 62,e > v « . (4.22) 

v m£n 

Once the voltage referred damage rate is obtained, the instantaneous damage rate can be obtained by simple 
application of the chain rule 

. v dD dD dv 

* (0a 1F= *a = *< v( '»*<<>- (4 ' 23) 
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We require of our model that damage, 23(f), be a monotonically increasing function, that is, the model will be 
conservative m that, healing of the battery will not be not allowed. Therefore we require the instantaneous 
damage rate to always be positive, thus equation (4.23) is modified to 

£(,> -f -f HH«'»i*«>i- 

Substituting the voltage referred damage rate into this equation gives the instantaneous damage rate as 

23(f) = 5222 x lO^c, e 4 - 621 ^^ |v(f)| . (4.25) 
Choosing v^l^V, gives at last 

23(f) = 5.222 x 10-^c, c «»^"> |y(f)| . ( 4. 26 ) 

^ This is the continuum damage model based on the Green-Hoffinan cyclic damage model. 

Q 

W Modified Damage Model 

U1 The Green-Hof&ian exponential model, equations (4.12) and (4.13), implies that cyclic damage 

U| occurs when the depth of discharge is zero. This is not physically plausible for an operational battery damage 

Ftf model. The degradation that occurs during shelf storage is not the same as cyclic damage. Thus we will 

*l modify the Green-Hoffinan model so that cyclic damage goes to zero as depth of discharge goes to zero. For 

* the model to allow zero damage at zero DOD we require the cyclic damage to have the form 
£? 

m 



(4.27) 



B . ^=c 2 (^>^)^c 2 (e^-l)=^ 

The constants c 2 and % are determined numerically to best fit the data presented in the table at the beginning 
of this section. This process yields the damage per cycle as 

= 1.0404 xlO" 5 ( e ^ D0D -l)=-L. (4.28) 

N f 

For DOD =35%, then this equation gives D^ c = 2.6316x1 0 s . To get the instantaneous damage, we will again 
assume a linear relationship between DOD (the charge amplitude) and cyclic amplitude of voltage deviation, 

DOD=c x v a . (4.29) 

To obtain c x , we adjust it until the simulated damage per cycle for a constant-plus-taper charge is equal to 
that given by Equation (4.28) for 35% DOD. Finding c x must wait until the end of the process however. 
This equation can be substituted back into the damage per cycle 

= 1.0404x 10- 5 (e 36Q2 ^ -l) (4.30) 
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Then, the integral of voltage referred damage rate, S(v) , for one charging cycle must be equal to the damage 
per cycle, and assuming damage formation only during the during the charging phase of the cycle 



p(v)dv = = 1.0404xl0- 5 ( e ijsa2 * v ° -l). (4 31) 
where 

= ~dv C4.32) 

is the voltage referred damage rate and Z>«) is the accumulated damage. Changing me limits of integration 
or Equation (4.31) gives 

J ^(v+v,^)^ = ]s(y+ Vmia )dv = 1.0404x10-* (e" 02 ™ -l) . (4.33) 
0 0 

The shifted voltage referred damage rate 5(y + ) is now determined to be 

<*(v + v^) = 1.0404 x 10~ s (3.602c,)e ,fi02c « v (4 34) 

or 

*(v) =3.7475x10^ Cl , v>Vaia , (4>35) 

This voltage referred damage rate exactly satisfies Equation (4.31). Having the voltage referred damage rate, 
the instantaneous damage rate can again be obtained by application of the chain rule 

m s ~d7 = * = *W«W0 • (4-36) 

Further requiring that the instantaneous damage rate is always positive gives 

*»-f-f|S|-*«w- «- 37 > 

Substituting the incremental damage rate into Ibis equation gives the instantaneous damage rate as 

£>(/) a 3.7475X10- 4 c, e 3 ^^ v ™)|^ _ (4 38) 
Choosing v mi)J =1.2V, gives 

S 3.7475 x 10"* c, e 3 - 602 ^ |v(/)J . (4.39) 
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Now, adjusting c, until Hie simulated damage per cycle for a constant-plus-taper charge is equal to that given 
by Equation (4.28) for 35% DOD, that isi?^ - 2.6316KH)- 5 , yields c, =0.9707 . 

It should be noted that the instantaneous damage rate could also be expressed as a function of the 
current and rate of change of voltage with charge by continued application of the chain rule, 



a,.. dD dD 



dv\ 


dq 


dq\ 


dt 



= *(v(0) 



dv 
dq 



m- (4.40) 



Here dv/dq is the slope of the charging curve on the v-q diagram. 

In deriving Equation (4.37), we required that the instantaneous damage rate always be positive. This 
was accomplished by using the absolute value of the voltage rate. This is justifiable in that whenever the 
voltage rate becomes large positive or negative, there are significant microscale changes taking place in the 
cell An alternative result would be to have the damage rate go to zero whenever the voltage rate goes 
negative. This is justifiable in that the cell may not be undergoing damage when it is relaxing, that is, when 
m the voltage rate is negative. Clearly, this topic requires further study. The damage rate in this case becomes 

H D(t) = 3.7475X10- 4 c x v(/), v(*)£0 

B ^) = o, v (0 <o (4 - 41) 

fjj In this case, the constant c x = 1.0901 . 

M- . 

p From the various forms for the damage, Equations (4.36>(4.41), it can be seen that changing voltage 

g£ is indicative of damage in the cell. The rate of change of voltage is an indication that the chemical reactions 

C?P are proceeding in a manner to cause damage on a microscopic level. Thus keeping the voltage from changing 

□ abruptly should allow longer life in these cells. This would seem to indicate that a "soft take-ofi" and a "soft 

£t landing" are important in the charging process to avoid damaging the cell. Looking back at the energy 

m balance equations with or without 02 production, Equations (4.10) or (4. 1 1) give the following expression for 

jfjj the voltage rate 



Mrt- 1 <*(r(QAS(Q) di(t) 

( 0 ~ Hf dt + IT ' (4,42) 

In words this equation says: 

changing terminal voltage oc changing heat of reaction + changing heating due to current, 

that is, the heating rate is intrinsically tied to the damage process. Stated differently, or on a microscopic 
scale, the equation indicates 



changing terminal voltage oc changing microscopic reordering rate + changing material flow rate. 

This is because the current is proportional to the concentration as indicated by equations (3.3) and (3.4). Both 
of these processes involve the rate of change of material, that is, material motion and acceleration, which 
would be related to internal damage rate. This would seem to infer that damage rate and internal loss rates are 
related, thus justifying the use of voltage-difference rate as the continuum damage rate model stress variable. 
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Now that we have an instantaneous damage rate model, we can use it to design an optimal charging 
profile that m ini m ises the accumulated damage per cycle. This will be the topic of Sections. 

4.6 Effect of Damage in the Essenttalized Model 

, ™ s secti <m will discuss the incorporation of die accumulated damage into the dynamic performance 
model First it is important to associate the accumulated damage with some physical mechanism. We will 
associate the accumulated damage with the inability to store charge, or to accumulate the stored material 
Thus this could be associated with the formation of the r - phase NiOOH, which actually does reduce the 
possible total stored charge. Likewise, overcharging always causes heating, which leads to the formation of 
/ - phase NiOOH. Consequently, we will associate the accumulated damage with a reduction in the ability to 
accumulate the stored material. 

The inability to store charge could potentially be manifested into all of the parts of the performance 
model; Faraday's Law, the electrode equation (or the linear-in-the-parameteis form), the charge storage term, 
the diffusion term, or the self-discharge. As we are associating the damage with a reduction in the inability to 
accumulate the stored material, we will focus on incorporating the damage into the charge storage term. 
Equation (3.22) becomes - - ° 

stored material with self-discharge and damage: c s (t) = i(t) —c s (t) - D(t) . 



p (4.43) 



But, as we have used the linear-in-the-parameters electrode equation, it must also be modified to indicate a 
loss of capacity. Thus Equation (3.24) becomes 

electrode equation with damage: v = ^ +k 2 ln(l + |i |)sgn(Q+£ 3 ln(c d )+ k 4 in((l-D(/))-e,), (4.44) 



5.0 Damage Mitigating Control 

5.1 Control Philosophy 

Control system design is usually separated into at least two phases. Phase one is the design of 
optimal trajectories and associated inputs, that move a given plant from one operating condition to another, 
while minimizing some performance measure. This subject is } die topic of the present section and the next 
Phase two is the design of a trajectory following controller (sometimes called a regulator or tracker) that 
provides a real-time control input perturbation to keep the plant operating near the designed optimal 
trajectory. This is discussed in Section 7. 

Our control philosophy is to charge the NiH2 cell in such a way that the damage incurred during the 
charging period is minimized, thus extending its cycle life. This is generally considered damage mitigating 
control or life-extending control. We will call the optimal charging profile life-extending charge. A 
damage mitigating control for use during discharge could also be considered, but the required discharge 
current was dependent on time-varying mission needs, and would require hardware modifications that were 
considered outside the scope of this research. Some discussion of possible damage mitigating strategies for 
discharge is included in the discussion of Section 8. 
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5*2 Performance Measure 



The performance measure to be minimized is the accumulated damage per recharge cycle. This can 
be obtained by integrating instantaneous damage rate with respect to time: 



(5.1) 



where D(i) is obtained from either of the damage models, Equations (4.26) or (4.39) , and t f for our problem 
is 54 Min. As it is, this performance measure can be minimized by using zero current, which does not 
recharge the cell. Thus, the performance measure must contain some information on the total charge required 
during charging. This is easily incorporated into a final state weighting term: 

J = w fi {c;(t,)-c,(t f )f + \t){t)dt, (5.2) 

0 

where c s (t f ) is the desired stored charge at the end of charge ( ^1), w fi is the cost weighting associated 
with the error in stored charge at the end of charge, and again D(t) is obtained fiom either of the damage 
models, Equations (4.26), (4.39), or (4.41), and //for our problem is 54 Min. Note that only one weighting 
term is necessary as it can be either greater or less than unity, but is constrained to be positive. In the 
implementation, we have chosen the self-discharge resistance to have a value which forces the final charge 
using constant-plus-taper charge (104% of 1420AMin) to have a value of 3900AMia The optimizer forces 
the life-extending charge to also have this final value. The constant charge (also 104% of 1420AMin) will 
have a final charge larger than this since the self-discharge current is slightly less than the constant-plus-taper. 

5.3 Original Optimization Approach 

To determine the optimal charging profile, or life-extending charge, that minimizes the performance 
measure, Equation 5.2, we will use a numerical optimization approach. The algorithm used is outlined in 
Figure 5. 1 . Some typical charging results are shown in Figures 5.2 a-d. 



Load previous best 
Setj=l 

Guess u(jT) 
Initialize simulation 

Main simulation loop 

Derivative function evaluations 
Numerical integration (Euler, T=l) 
Damage rate calculation 
Integrate damage rate 
Bid main loop 

Include final state error into Performance Measure 
Change u(jT) using gradient search and repeat until optimal 
Increment j and repeat 
Reset j=l and repeat until optimal 
Store present best 



Figure 5.1. Outline of original optimization approach. 
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5.4 Charging Profile Design for Specific Battery 

Hie optimization originally was accomplished by sampling the charging current in time, and 
optimizing each little sample increment This optimisation approach is veiy slow and requires many hours of 
computation time. The results shown in the last section, however, seemed to indicate that the optimal 
charging profile was similar to a parabola. Using this idea allows a much fester optimization to be 
accomplished, although the results are probably a little conservative as the truly optimal currents at each point 
in time are not being calculated. The fast optimization approach is outlined in Figure 5.3. Here a ten point 
cubic spline is chosen to represent the current profile so as to allow some variation from a parabola (other 
time functions were also considered which gave similar results). Charging profiles are shown in the next 
several plots. In performing the optimization, it has been important to be sure to include the cost of requiring 
the zeroing of a non-zero current at the end of charging, as the large voltage rate there may incur significant 
damage. Also notice the constant current charge is predicted to last for only 25465 cycles, while the actual 
cell is still cycling after 50000 cycles. We believe this indicates either some problem with the cycle life data, 
or that this particular cell has used a design different from that of the cycle life data. 

Recognizing that the cell voltage response is diffusion limited on discharge, along with the 
knowledge that much of the damage occurs near the end of charge when most of the overcharge occurs, 
allows us to try a charge-discharge profile that backs off from the 100% charged point Said otherwise, we 
are going to determine the life-extending charge profile from 60% to 95% fully charged. On discharge, the 
cell will still use 35% of the AHr capacity since it is diffusion limited at cutoff. On charge, the cell will be 
kept away from the overcharging region, and thus reduce the damage per cycle. The resulting life-extending 
control will be implementable in real-time due to the use of an on-board dynamic observer which can keep 
track of the stored charge material. The use of a real-time observer is discussed in the next section. 

The expected improvement associated with the life-extending charge is summarized in Tables 5.1, 
5.2, and 5.3 which show a significant potential extension in life expectancy. The percentage life increase 
should be regarded as approximate due to the use of the fast optimizer and relatively large simulation 
timestep. It is interesting that the optimal solution using the | dv I dt\ term in the damage rate model (Table 5.1 
and Figure 5.4), gives the same optimal charging curve as incurring damage only when dv/dt is positive 
(Table 5.3 and Figure 5.6), although the actual cycle life is different because of different damage rate models. 



Load previous best 

Guess u(t) function coefficients 
Initialize simulation 
Main simulation loop 

Derivative function evaluations 
Numerical integration (Euler, T=l) 
Damage rate calculation 
Integrate damage rate 
End main loop 

Include final state error into Performance Measure 
Change u(t) function coefficients using gradient search and repeat until optimal 
Store present best 



Figure 5.3 Outline of fast optimization approach. 
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Life-extending-Charge 


ox t HS» 

/o i*iie 
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Damage per 
65-100% Cycle 


Cycles to 
Failure 




Constant-Current 
.Charge 


3.9269e-005 


25465 


2.0780e-005 


48123 


88.98% 


Constant -r laper 
Charge 


2.6316&-005 


38000 


2.0780e-005 


48123 


26.64% 


Table 5.1 


Life-extending charge compared to constant current and constant current + taper. 




Damage 
per Cycle 


Cycles to 
Failure 


Life-extending-Charge 


%Life 
Extension 








Damage per 
60-95% Cycle 


Cycles to 
Failure 




Constant-Current 
Charge 


3.9269e-005 


25465 


1.9535e-005 


51190 


101.02% 


Constant + Taper 
Charge 


2.6316e-005 


38000 


1.9535e-005 


51190 


34.71% 


Table 5.2 60%-95% Life-extending charge compared to constant current and constant current + taper. 




Damage 
per Cycle 


Cycles to 
Failure 


Life-extending-Charge 


%Life 
Extension 








Damage per 
65-100% Cycle 


Cycles to 
Failure 




Constant-Current 
Charge 


3.9269e-005 


25465 


2J2080e-005 


| 45290 


77.85% 


Constant + Taper 
Charge 


2.6316e-005 


38000 


2.2080e-005 


45290 


19.18% 



Table 5.3 Life-extending charge compared to constant current and constant current + taper using only 
positive dv/dt in the damage rate law.. 
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Figures 5.5 a&b. Comparison of 60%-95% and standard charging methods, charge and current 
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Figures 5.6 c&d Comparison of life extending (+dv/dt damage only) and standard charging methods. 
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Figures 5.6 e&f. Comparison of life extending (+dv/dt damage only) and standard charging methods. 
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6.0 Suggested Real-time Implementation 

6.1 Real-time Observers 



no71Y > ^ m Portant result of state space control theory is the concept of a state observer (Luenberger 
(1971)). The total number of dynamic "states," or "state variables," in a system is equal to the number of 
energy storage events in the system, which is also equivalent to the number of time derivatives that are 
I S«S I ^^TSS e S? ions ^Presenting that system (Freidland(1986)). The "state" of a given system 
is the collection of aU the "state variables" of that system at a given point in time. This coUection is all the 
mfonnation needed to advance the "state variables" forward in time through the differential equations 
representing the system. With these definitions, the concept of a state observer will now be clarified 

11 i A -f!? obser Y er is effectivel y a dynamic model of a given system that runs in real-time, and in 
parallel with the actual system. Typically, only a small number of the dynamic states of an actual system are 
measureable, or are sensed. In the case of an electrochemical cell, typically only terminal voltages and 
currents, and perhaps some pressures and temperatures are all that are sensed. The state observer is then used 
to provide an estimate of the non-sensed variables. There are usually two forms for state observers, the open- 
ly loop observer and the closed-loop observer. 

S ^ open-loop observer receives the same inputs as the actual system, and then provides the collection 

m of dynamic states as an output (Jacquot(1981)). If the model being used in the observer is accurate, then all 
yi toe dynamic states from the observer will accurately represent the true states in the actual system. This has 
m been applied successfully to a jet engine (Merrill and Lorenzo (1988)), as well as other places. A closed-loop 
M observer receives the same inputs as the actual system, as well as all of the sensed dynamic states and outputs 
M the actual system (Jacquot (1981)). Comparing these sensed variables with the corresponding simulated 

a variables, the difference can be used as a feedback signal to force the observer dynamic states to converge to 
p the true dynamic states of the actual system (see Figure 2). This approach does not require the model to be as 
gjf accurate as m the open-loop observer due the correction of the feedback from the sensed variables. Closed- 
0 Io °P observers are used in most state space control designs (Maciejowski(1989), Stengel (1986)), where they 
^g*. have also been used to filter noise out of the measurements. 

0 

fP A real-time observer for an electrochemical cell incorporates the hybrid cell model, as discussed in 

the previous section, into a microprocessor based simulation that would be required to run in real-time. 
Sensed variables from an actual cell would then be input to the microprocessor simulation to force the 
dynamic states of the real-time cell observer (simulation) to converge to the corresponding states in the actual 
cell. The feedback is crucial for this problem as the actual system is highly nonlinear as well as spatially 
distributed, and the model dynamic states represent aggregations of many actual states. Even so, the feedback 
to the observer will force the errors in the measured variables to zero, and the unmeasured variables will 
accurately represent the corresponding actual variables (Krob (1992)). It is expected that a real-time observer 
for an electrochemical cell will be able to provide accurate information about state of charge, active species 
concentrations, and internal temperature. Any theory required for the fractional order components in the 
observer will need to be developed. More advanced observer designs may also provide some indication of 
damage as discussed below. 
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Figure 2. Real-time observer structure. 



6.2 Observer Based Control 

T^picaUy real-time observers are used in two ways, for continuous-state feedback and for decision 
malomg. Much of modern control theory is based on state feedback. State feedback requires knowledge of all 
the dynamic states of the given system, which are then fed back to the system inputs to improve overall 
system performance. As most of these dynamic stales are not measurable, most state feedback controllers use 
the estimated states from real-time observers for feedback purposes. 

Observers can also be used as part of the general decision making process. To do so, then requires an 
extra layer of intelligence. For example, if the observer predicts that certain dynamic states have varied too 
fir from where they are supposed to be, the higher level of intelligence can determine if some sensors have 
felled, if some actuators have foiled, or if the actual system has become damaged in some way (Merrill and 
Lorenzo (1988)). Of course, this intelligent layer must be programmed by a knowledgeable user. Another 
decision making application could be to shut the system down if certain observed dynamic states exceeded 
those which would normally be deemed healthy for the actual system. These decision making control 
strategies will be addressed later when discussing advanced controllers and life-extending control. 

Some of the issues associated with observer based control of batteries and fuel cells require a* clear 
understanding of improved performance. Some potential definitions of improved performance might include 
maximum energy out, maximum power out, constant voltage out, constant current out, constant power out, 
prolonging state of health, or minimum damage rate, among others. Other issues are more hardware related. 
For example, in a collection of cells, control of each cell would probably be important An array of observers 
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one for each cell, would then be able to provide information about individual cells. As long as each cell could 
be independently controlled, this information could actually be used handle any degradation or failures of 
mdmdual cells. Of course this also requires the usage of power electronics to control the charging and 
^charging of individual cells (Tarter (1993)). This could include the use of buck-boost DC-DC carters 
to provide constant voltages or currents out, as cell voltages and states of charge change. Finally, any theory 
necessary for control of the ftactional order components will be developed ~ 

63 Advanced Real-time Observers 

An advanced real-time observer for electrochemical cells would require an extra level of intelligence 
over the standard real-time observer. Specifically, the algorithm used for the identification of the cell hybrid- 
model parameters would need to be automated (Goodwin and Sin(1984)). This would then allow two things. 
First, it would allow the tracking of parameters as a cell changes during its lifetime. Some knowledge of 
parameter range could indicate when a cell was nearing failure. Second, it would allow the on-line 
determination of the cell model parameters to begin with, or effectively, the correction of the initial 
parameters provided to the real-time observer. Effectively, the real-time observer could then "learn" any 
battery. An advanced real-time observer should also include models of any damage mechanisms as a function 
of the system states. 

It is felt that degradation algorithms should be included into an observer structure to give a real-time 
damage indicator for a given cell (Tobias and Trindale (1995)). 



6.4 Advanced Controllers 

An advanced controller for an electrochemical cell would also require an extra level of intelligence. 
Specifically, die advanced controller would use the information provided by the advanced real-time observer 
to maintain cell performance as the cell's parameters change with time. This would require a continuous 
change of controller gains as the effect of physical changes in a given cell were updated by the advanced real- 
time observer (Goodwin and Sin(1984)). 

Secondly, an advanced controller would need to implement some of the decision making skills as 
discussed earlier in the Observer Based Control section. If the estimates provided by the real-time observer 
vary too much from the actual cell outputs, then decisions would be required concerning sensor, actuator, or 
cell health. If any failures were detected, additional intelligence would be required to accommodate a 
detected failure. Of course, the chosen battery and control structure must allow accommodation. 

Finally, assuming some form of a damage model was available, the controller could then make 
performance-lifetime trade-off decisions. Associated with this intelligence, there may need to be off-line 
optimization studies performed that determine the most efficient charging (or discharging) profile which will 
also prolong die cell lifetime. This is related to the area of optimal control. One of the most well known off- 
line optimal control results is the minimum-time-to-climb jet fighter intercept solution which was not at all 
expected (Bryson and Ho(1975)). Another example of off-line optimization results is the reusable rocket 
engine damage mitigating control mentioned earlier (Lorenzo (1992)). Ultimately, designing a cell charging 
(or discharging) profile that extends the useful life of the cell by minimizing accumulated damage would be a 
useful result. 
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Figure 3. Advanced control system using advanced real-time observer. 
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7.0 Summary and Future Research 

7.1 Summary 

^ In this report, we have presented a methodology for extending the life of N£H2 cells via a damage- 
mitigatog charging current This has required the development of an essentialized dynamic model that will 
accurately reconstruct internal states in a given cell The parameters of this model are readily obtainable from 
either off-line data, or adaptively on-line from real-time data. To determine the life-extending charging 
current, an instantaneous damage model that determines the damage rate at any point in time based on voltage 
measurements, has been developed. Using the instantaneous damage rate model with the essentialized cell 
model, allows the development of a the damage mitigating charge profile for extending the life of the cell. 
This charge profile is obtained via optimization, and can be modified by trading recharge for damage. Once 
the life^xtending charge profile has been obtained, the implementation of this profile is performed via a real- 
time controller that uses a real-time observer for reconstructing the internal cell variables based on current and 
voltage measurements at the cell terminals. This real-time controller implementation can be used over the 
entire life of the cell by updating any combination of the model parameters, control parameters, or charging 
profile based on the data from the present charge-discharge cycle. ° 

CJ Model Development 

jg 1 ■ Developed an essentialized NiH2 cell model structure. 

fS 2 - Developed a method for determine the essentialized model parameters for charge- 

discharge data. 

ffi 3. Developed a deeper understanding of the charge-discharge hysteresis behavior in the 

y voltage-charge plane. 

j*f Damage Model 

g? 1 ■ Developed an instantaneous damage model for NiH2 cells based on cell voltage 

f measurements. 
2. The instantaneous damage rate is shown to depend on voltage and voltage rate. 
n 3. It is shown that the size of the voltage-charge hysteresis curve is related to the damage 

{Hi per cycle. 

Life-extending Charge Development 

1 . A life-extending charging profile has been developed by optimizing the charging current 
with respect to the damage per cycle and total recharge. 

2. This life-extending current should provide significant life extension over the present 
recharging methods. 

3. Significantly more life extension for a given cell can be obtained by not recharging the 
cell to 100%. 

4. Some heuristic rules for minimum damage charging that have emerged are: 

a. Increase the charging current slowly at the beginning of charge. 

b. Put more of the charge in during the first half of charging. 

c. Decrease the current to zero slowly at the end of charge. 

d. During charge, try not to allow the terminal voltage to increase much beyond the 
thermoneutral voltage. 

e. Try to keep the voltage-current hysteresis-curve area as small as possible. 
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Observer Development 



L ^vdSed 1 " 0 ^ 61 ° bSerVer theory ( necessai y f ° r control implementation of the model) has been 

2 - ^parameter identification technique is discussed for allowing a real-time essentialed model to 
xracK changes in its parameters in real-time over the life of a celL 

l** 301 ? 6 C<mtr0ller for implementing the optimal charging profile in real-time based on 
observed system states is discussed. on 

4. A learning controller that continuously updates the charging profile based on changing 
cell parameters has been discussed 

7.2 Future Diections 

Overcharge Control 

H 2M Ju£?n!ft!?^ By monitoring the 

2e P r + , ^ a 1 MH2 ^ rt be P° ssible to provide some control for overcharging. The cell pressure 
was not mcluded m our essentialized model, but could be without much increase in complexity. We could 
then relate then H2 pressure to both overcharge and damage. P * 

By using our instantaneous damage model, both a real-time damage rate sensor and an on-line 
damage monitor could be developed The damage rate sensor would require knowledge of the amount of 
stored material, and the present voltage. The damage monitor, or life predictor, could be developed using the 
damage rate sensor, or by using an accelerated degradation model as discussed in Tobias and Trindade. 

. Ultimately a charge-discharge controller could use feedback from these sensors to trade-off lifetime 
with performance. 



Other Chemistries and Fuel Cells 

It is expected that this methodology can be readily applied to other chemistries, such as Li-based cells 
if cycle-life data is available. If this data is not available, a cycle-life model can be inferred from an 
accelerated degradation approach (Tobias and Trindade). 

With regard to fuel cells, mis methodology should also be applicable, however it will be necessary to 
modify the essentialized model structure for fee two feel sources, ft will also be necessary to understand the 
damap mechanisms associated with fuel cells. Alternatively, it is expected that a tracking controller should 
be able to immediately improve performance in a fuel cell. 



Discharge Damage Control 

Although discharge control was considered outside the scope of this research, various possibilities exist to 
reduce damage during the discharge part of fee cycle. It was assumed in fee damage modeling that the entire 
damage during fee cycle took place over the charging part of the cycle. It is of course possible to apportion 
feat damage over both the charge in the discharge parts of the cycle. In the former, fee damage model is 
dictated by fee form of the damage per cycle and would remain very much the same. Future research could 
generalize the damage model to include damage on both charging and discharging With fee availability of 
the generalized damage model, various strategies could be investigated which would control damage during 
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discharge The additional aspect of such research is that it could influence the manner in which the battery 
power is allocated. The acceptability of this would depend on the specific application. 

Damage Testing and Future Data Needs 

The present research allowed access to the Crane database. From our observation of the database, which is 
certainty not a detailed analysis of the database, we have found it to be an excellent resource for evaluation of 
battery designs and battery chemistries, etc. However, it was not clear mat appropriate data was available in 
the database, for evaluating and creating damage models for particular battery types. It has been our 
observation, that in some of the battery literature, me battery design aspects and the battery operational 
aspects are not clearly delineated. It is important to separate these two areas, and to develop as a separate 
discipline, battery operational theory. In particular, clean operational data seems to be required. By mis we 
mean tests to failure, with key parameters held constant One direction of future research would be to 
establish testing methods specific to obtaining usable damage (operational) models. This involves accurate 
assessment of the equation forms for cycle damage as well as proper statistical treatment of the data. 
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CLAIMS FOR UA.439P ENTITLED OPTIMAL BATTERY CHARGING 
FOR DAMAGE MITIGATION 

What is claimed is: 

1. A method for charging an electrical storage device so as to extend the life thereof, 
comprising: 



charge plane. 

2. The method according to claim 1, further comprising: 

measuring voltage values of said structure based upon said charge-discharge 
hysteresis; and 



3. The method according to claim 2, further comprising: 

developing a charging profile based upon said instantaneous damage rate, 
wherein said charging profile optimizes a charging current with respect to the damage 
per cycle so as to "extend the overall life of the electrical storage device. 



developing an essentialized cell model structure of the electrical storage device; 
determining model parameters for charge-discharge data of said structure; and 
determining charge-discharge hysteresis behavior of said structure in a voltage- 



deriving an instantaneous damage rate from said measured voltage values. 
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